Pion productions by proton and Helium-3 on 197 Au target 
at beam energies of 2.8, 5, 10 and 16.587 GeV/nucleon 
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Based on a Relativistic Boltzmann-Uehling- Uhlenbeck transport model, proton and 3 He induced 
reactions on 197 Au target at beam energies of 2.8, 5, 10 and 16.587 GeV/nucleon are studied. It is 
found that compared with proton induced reactions, 3 He induced reactions give larger cross sections 
of pion production, about 5 times those of the proton induced reactions. And more importantly, pion 
production from 3 He induced reaction is more inclined to low-angle emission. Neutrino production 
via positively charged pion is also discussed accordingly. 
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There exists considerable interest in the possibility 
that one type of neutrino may transform into another 
type while propagation [ij. It was also argued that neu- 
trino oscillations are related to stellar collapse [2| and 
spontaneous parity non-conservation Q. Nowadays at- 
mospheric , reactor Q , solar neutrino @, 0] and ac- 
celerator neutrinos [Toj | provide compelling evidences for 
neutrino mass and oscillation. It is interesting to note 
that the precise neutrino oscillation parameters arc deter- 
mined by the KamLAND [ll| recently. For more about 
neutrino physics, please see Ref. 12]. In the accelerator 
based neutrino experiments, a key issue toward the de- 
velopment of a muon collider or neutrino beam based on 
a muon storage ring is the design of a target/capture sys- 
tem capable of capturing a large number of pions. These 
pions then proceed into a decay channel where the re- 
sultant muon decay products are harvested before being 
conducted into a cooling channel and then sub seq uently 
accelerated to the final energy of the facility [13]. Un- 
derstanding of the production of pions in proton interac- 
tions with nuclear targets is thus essential for determin- 
ing the flux of neutrinos in accelerator based neutrino 
experiments [HI, EH . A large amount of data were col- 
lected by the HARP Collaboration experiments for the 
above physical subjects recently [l6|. In fact, in acceler- 
ator based neutrino experiments, some times one needs 
Helium-3 induced reaction. There are many simulation 
methods focus on such studies, such as the phenomeno- 
logical Monte Carlo generators GEANT4 [l7j and MARS 
[18j and other theoretical works that considering more 
physical processes fillip . In this article, after check- 
ing the reliability a relativistic transport model (ART) 
we made comparative studies of pion production in pro- 
ton and 3 He induced reactions on 197 Au target at inci- 
dent beam energies of 2.8, 5, 10 and 16.587 GeV/nucleon. 
And finally, we simply discussed neutrino production via 
positively charged pion. 

The well-known Boltzmann-Uehling-Uhlcnbeck (BUU) 
model [2f| has been very successful in studying heavy- 
ion collisions at intermediate ener gies . The ART 
(A Relativistic Transport) model [27| is the rela- 
tivistic form of the BUU model, in order to simu- 
late heavy-ion collisions at higher energies, some new 



physics were added. It includes baryon-baryon, baryon- 
meson, and meson-meson elastic and inelastic scatter- 
ings. The ART model includes the following baryons 
N, A(1232), 2V*(1440), iV*(1535), A, S, and mesons 
7r, p, u), rj, K, as well as their explicit isospin degrees 
of freedom. Both elastic and inelastic collisions among 
most of these particles are included. For baryon-baryon 
scatterings, the ART model includes the following inelas- 
tic channels: NN O N(AN*), NN O A(AAT*(1440)), 
NN O NN(npuj), (NA)A O NN*, and AiV*(1440) O 
AW* (1535). In the above, N* denotes either AT*(1440) 
or AT* (1535), and the symbol (AN*) denotes a A or an 
N* . For meson-baryon scatterings, the ART model in- 
cludes the following reaction channels for the formation 
and decay of resonances: ttN O (AN* (1440) AT* (1535)), 
and rjN O AT* (1535). There are also elastic scatter- 
ings such as (ftp) (NAN*) -> (np)(NAN*). For meson- 
meson interactions, the ART model includes both elas- 
tic and inelastic 7T7t interactions, with the elastic cross 
section consisting of p meson formation and the remain- 
ing part treated as elastic scattering. Also included are 
reaction channels relevant to kaon production. The ex- 
tended ART model is one part of AMPT (A Multi-Phase 
Transport Model) model [28[. We use the Skyrme-type 
parametrization for the mean field which reads (27| 



U(p)=A(p/p ) + B(p/p y 



(1) 



Where a = 7/6, A = -0.356 MeV is attractive and B = 
0.303 MeV is repulsive. With these choices, the ground- 
state compressibility coefficient of nuclear matter K=201 
MeV. More details of the model can be found in the orig- 
inal Ref. 



27] 



To check the reliability of using the ART model to 
study the cross sections of pion production in proton or 
3 He induced reactions, we first made a comparison of 
pion production in p+Au reaction at an incident beam 
momentum of 17.5 GeV/c between the theoretical simu- 
lations and the E910 data [l4| as shown in Fig.Q] The top 
panel of Fig. [1] shows the inclusive differential cross sec- 
tions of pion production from p+Au at an incident beam 
momentum of 17.5 GeV/c. We can see that for both ir~ 
and 7r + , our results fit the E910 data very well, especially 
at higher momenta. Pion production of p+Cu reaction 
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FIG. 1: (Color online) Top: Production cross sections of 
7r~(left column) and ir + (right column) from p+Au at the 
incident beam momentum of 17.5 GeV/c (Et ea m ~ 16.587 
GeV/nucleon) shown in bins of cos# (relative to beam direc- 
tion). Numbers in the legend refer to the center of each bin. 
Data are taken from Ref. Bottom: Same as p+Au case, 
but for 3 He+Au. 



at incident beam momenta of 12.3 and 17.5 GeV/c also 
fit the E910 data [lH very well. From Fig. [TJ we can also 
see that the cross sections at low-angle (0.9 < cosO < 1) 
are evidently larger than those at high-angles, especially 
for energetic pion mesons. As a comparison, we also give 
the case of 3 He+Au at the incident beam momentum 
of 17.5 GeV/c as shown in the bottom panel of Fig. [T] 
From the bottom panel of Fig. [1] it is seen that differen- 
tial cross sections of pion production of the 3 He induced 
reaction are about 5 times those of the proton induced 
reaction at the incident beam momentum of 17.5 GeV/c. 
The cross sections at low-angle (0.9 < cosO < 1) are also 
much larger than those at high-angles. 

To make comparisons systematically between p+Au 
and 3 He+Au at different incident beam energies, we plot 
Fig. rjHH From the top panels of Fig. [50 we can clearly 
see that as incident beam energy decreases, cross sec- 
tions of pion production of p+Au also decrease. This is 
understandable since pion production mainly comes from 
decays of resonances and energetic nucleon-nucleon col- 
lisions give more resonances. We can also see that as 
beam energy decreases, the energetic pion mesons also 
decrease rapidly. The bottom panels of Fig. [2j4] are 
the cases of 3 He+Au at different incident beam energies. 
Also it is seen that as beam energy decreases, differen- 
tial cross sections of pion production decrease rapidly, 
especially for energetic pion mesons. Compared with 
p+Au, as the increase of incident beam energy, cross 
sections of pion production at low-angle of 3 He+Au in- 
crease more rapidly, especially for energetic pion mesons. 



FIG. 2: (Color online) Top: Production cross sections of 
7r~(left column) and tt + (right column) from p+Au at the in- 
cident beam energy of 10 GeV/nucleon shown in bins of cos 9. 
Numbers in the legend refer to the center of each bin. Bottom: 
Same as p+Au case, but for 3 He+Au. 




0.0 0.3 0.6 0.9 1.2 0.3 0.6 0.9 1.2 

p [GeV/c] 







5 GeV 

T 3 He +Au->«" + X 1 


Z 0.95 ' - . ~ 

- 0.85 '--v 

_ 0.75 

0.65 

_ 0.55 " _ 

: 3 He + Au > ti + + X 



O.O 0.3 0.6 0.9 1.2 0.3 0.6 0.9 1.2 

p [GeV/c] 



FIG. 3: (Color online) Top: Production cross sections of 
7r~(left column) and ix + (right column) from p+Au at the in- 
cident beam energy of 5 GeV/nucleon shown in bins of cos 9. 
Numbers in the legend refer to the center of each bin. Bottom: 
Same as p+Au case, but for 3 He+Au. 



In the incident beam energy region from 2.8 to 16.587 
GeV/nucleon, cross sections of pion production at low- 
angle (0.9 < cosO < 1) and high momentum of 3 He+Au 
are 5^10 times those of p+Au case. Using the AMPT 
model we also made simulations for p+Au at incident 
beam energies from 50 to 100 GeV/nucleon, cross sec- 
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FIG. 4: (Color online) Top: Production cross sections of 
7r~(left column) and ir + (right column) from p+Au at the 
incident beam energy of 2.8 GeV/nucleon shown in bins of 
cos 9. Numbers in the legend refer to the center of each bin. 
Bottom: Same as p+Au case, but for He+Au. 
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FIG. 5: (Color online) Angle distributions of pion multiplicity 
of p+Au and 3 He+Au at different incident beam energies, in 
bins of cos 6. 



tions of pion production at low-angle (0.9 < cosQ < 1) 
are both about 20 times those of p+Au at incident beam 
energy of 16.587 GeV/nucleon, indicating the saturation 
of cross section of pion production at the beam energy of 
about 50 GeV/nucleon. 

We next turn to the study of angle distributions of pion 
multiplicity of p+Au and 3 He+Au at different incident 
beam energies. From Fig. [5j we can see that whether 
for p+Au or 3 He+Au, pion emission at low-angle in- 
creases rapidly, especially at higher incident beam en- 
ergies. From these plots, we can also see that the low- 
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FIG. 6: (Color online) Angle distributions of charged pion 
relative multiplicity (N coa g/N tota i) for p+Au and 3 He+Au 
reactions at different incident beam energies, in bins of cos 9. 



angle's pion emission is more pronounced for 3 He+Au 
than p+Au, especially at higher incident beam energies. 
Pion numbers from 3 He+Au at low-angle (0.9 < cosO < 
1) are about 5 times those of p+Au. This indicates that 
3 He+Au at high incident beam energy is more suitable 
for neutrino experiments compared with p+Au. Fig. [S] 
shows angle distributions of pion relative emitting num- 
ber at different incident beam energies. It is seen that 
at the high incident beam momentum 17.5 GeV/c, rel- 
ative emitting number at low-angle (0.9 < cos8 < 1) 
can reach about 50% for 3 He+Au. While at incident 
beam energy 2.8 GeV/nucleon, relative emitting number 
at low-angle (0.9 < cosO < 1) reaches only about 25%. 
At the studied beam energy region, we can clearly see 
that the 3 He induced reaction on Au target causes larger 
proportional low-angle pion emission, especially at higher 
incident beam energies, about 5%^ 10% larger than that 
of the proton induced reaction on Au target. 

Fig. [7] shows cross sections of charged pion produc- 
tion of p+Au and 3 He+Au reactions at different incident 
beam energies. We can see that at the incident beam 
energies studied here, cross sections of pion production 
of 3 He and proton-induced reactions on target Au in- 
crease about 3 times. The cross sections of 3 He-induced 
reaction on target Au at the incident beam energy of 2.8 
GeV/nucleon are larger than those of p+Au reaction at 
the beam energy of 16.587 GeV/nucleon. Fig. [5] shows 
the ratio of cross sections of charged pion production of 
p+Au and 3 He+Au reactions at different incident beam 
energies. We can clearly see that cross sections of charged 
pion production from 3 He induced reaction are about 5 
times (larger than A3h £ /Ah = 3) those of the proton in- 
duced reaction. 

Since the work focuses on the proton versus 3 He re- 
sults, one wonders what scaling behavior with projectile 
nucleon number would one expect from a "standard" cas- 
cade model (without mean-field modifications)? Fig. [5] 
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FIG. 7: (Color online) Cross sections of charged pion pro- 
duction of p+Au and 3 He+Au reactions at different incident 
beam energies. 



FIG. 9: (Color online) Ratio of mean pion production per nu- 
cleon of projectile with mass number A and 1 at the incident 
beam energy of 2.8 GeV/nucleon. The target is 197 Au and the 
impact parameters are set to be and random, respectively. 
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FIG. 8: (Color online) Ratio of cross sections of charged pion 
production of p+Au and 3 He+Au reactions at different inci- 
dent beam energies. 



shows the ratio of mean pion production per nucleon of 
projectile with mass number A and 1 (proton) at the 
incident beam energy of 2.8 GeV/nucleon (by the ART 
cascade model). We can see that pion production per 
projectile nucleon is roughly the same with different pro- 
jectile mass number A, i.e., the produced total pion num- 
ber is roughly proportional to projectile mass number 
A. This indicates each nucleon in the projectile excites 
pion production almost dependently. But for the projec- 
tile which mass number is smaller than the target mass 
number, the scaling behavior that total pion number is 
roughly proportional to projectile mass number A is not 
strictly correct. In fact, the ratio of — /A is always larger 
than 1, as shown in Fig. [HI This is because each nucleon 



in the projectile does not excite pion production indepen- 
dently. One nucleon in the projectile may give energy to 
the nucleon in the target, but does not produce pion. The 
other nucleon in projectile may also collide with the nu- 
cleon which has obtained energy from the other nucleon 
of the projectile. Thus the probability of pion produc- 
tion accordingly increases for the other induced nucleon 
in the projectile. This correlation of different incident 
nucleons of the projectile does not increase linearly with 
projectile's mass number due to the marginal collision of 
the induced nucleon. Thus we see about 5 times pion 
production of the 3 He induced reaction compared with 
the proton induced reaction. For random impact param- 
eter case, smaller /A is due to the marginal collisions 
of the induced nucleon in the projectile. 

To demonstrate neutrino production by proton or 3 He 
induced reactions, we plot Fig. [TUJ angle distributions of 
neutrino via ir + decay in proton and 3 He induced reac- 
tions on target Au. Assuming 7r + — > fi + + and pion 
decays into \i and isotropically in its frame of refer- 
ence, we can thus obtain neutrino distribution by assum- 
ing it rest mass 1 eV. Fig. [TUl shows angle distributions of 
neutrino production from positively charged pion dacay 
in p+Au and 3 He+Au reactions at incident beam ener- 
gies of 2.8 and 10 GeV/nucleon. We can clearly see that 
neutrinos from 3 He induced reaction are more inclined 
to low-angle emission than proton induced reaction, this 
situation is clearer for high incident beam energy. Be- 
cause the energy distribution of the emitting neutrinos 
are important for neutrino-nucleus experiments [29|-[32j], 
we also plot the energy distributions of the produced neu- 
trinos at low and high angles as shown in Fig. QTJ We 
can see that the produced neutrinos possess different en- 
ergies from about 1 MeV to 1000 MeV and more. The 
most probable energy is about 30^70 MeV for several 
GeV incident beam energy. Moreover, we can see that 
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FIG. 10: (Color online) Angle distributions of neutrino pro- 
duction from positively charged pion decay in p+Au and 
3 He+Au reactions at incident beam energies of 2.8 and 10 
GeV/nucleon, in bins of cos#. 
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FIG. 11: (Color online) Energy distributions of neutrino pro- 
duction at angles 0.9 < cos9 < 1 and 0.5 < cos6 < 0.6 from 
positively charged pion decay in p+Au and 3 He+Au reactions 
at incident beam energies of 2.8 and 10 GeV/nucleon. 



neutrinos from low-angle possess more energy than those 
from high-angles. Note here that neutrino production 
can also from other channels [33[ , especially for energetic 
collisions. For physical experiments relevant to neutri- 
nos, detailed studies of the numbers, the energy spectra 
as well as the species of emitted neutrinos are very nec- 
essary and therefore the simulations related to neutrino 
production are also become important [33[. 



In conclusion, proton and 3 He induced reactions on 
197 Au target at beam energies of 2.8, 5, 10 and 16.587 
GeV/nucleon are studied in the framework of the Rel- 
ativistic BUU transport model. It is found that com- 
pared with proton induced reactions, 3 He induced reac- 
tions give larger cross sections of pion production, about 
5 times those of the proton induced reactions. And 3 He 
induced reactions are more inclined to low-angle's pion 
emission. Simulations demonstrate that neutrino emis- 
sion via positively charged pion decay is also inclined to 
low-angle emission. 
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